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Abstract 

> . 

We establish anomalous inertial range scaling of structure functions for a model of 
advection of a passive scalar by a random velocity field. The velocity statistics is 
\q . taken gaussian with decorrelation in time and velocity differences scaling as |x| K//2 in 

' space, with < k < 2. The scalar is driven by a gaussian forcing acting on spatial 

scale L and decorrelated in time. The structure functions for the scalar are well 
defined as the diffusivity is taken to zero and acquire anomalous scaling behavior 
for large pumping scales L. The anomalous exponent is calculated explicitly for the 
4 th structure function and for small k and it differs from previous predictions. For 
all but the second structure functions the anomalous exponents are nonvanishing. 
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In 1941 A.N. Kolmogorov argued that in fully developed turbulence a range of spatial 
scales exists where the velocity structure functions acquire a form independent of the IR 
and UV cutoffs provided by the scale of energy pumping and dissipation respectively. 
Ever since a debate has been going on as to whether there are corrections to the scaling 
exponents predicted by Kolmogorov and whether such corrections depend on the dissipa- 
tion or the pumping scale or both. This question being still quite open for Navier-Stokes 
turbulence both experimentally and theoretically, it is useful to consider it in the context 
of simpler models that are nevertheless expected to display phenomena similar to the 
Navier-Stokes equations. 

In this letter we consider one such model that has attracted much attention recently 
[0, fj, [5], f|, []], namely that of passive advection in a random velocity field v(£, x) of a 
scalar quantity T whose density T(t, x) satisfies the equation 

d t T+(v-V)T-vAT = f, (1) 

where v denotes the molecular diffusivity of the scalar T and f(t, x) is an external source 
driving the system. 

We take v(t, x) and /(£, x) to be mutually independent Gaussian random fields with 
zero mean and covariances 

<^(t,x)^'(t / ,x') >= 8(t - if) D ij '(x - x') , (2) 
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< f(t, x) f(t', x') >=8{t- C(^) = - Cl(x - x') 



(3) 



Both are thus decorrelated in time, a fact that leads to an exact solution for the correlation 
functions of the scalar. The forcing covariance C is assumed to be a real, smooth, positive- 
definite function with rapid decay at spatial infinity so that the forcing takes place on the 
("integral") scale L. 

The velocity covariance D is taken to mimic the situation in real turbulent flow with 
structure function < (v(i, x) — v(i,0)) 2 > proportional to |x| K for k > 0. Concretely, we 
set 

D»(x) = D J e ik x (k 2 + m 2 )- (3+K)/2 (5 ij - W/k 2 ) ^ , (4) 

where the transverse projector in the Fourier space ensures the incompressibility of v. 
Small m 2 is an infrared cutoff making the integral convergent for < k < 2. Writing 
£>(x) = D(0) - Z>(x) , we have 

= uJrfLyv ^D m-, (5) 
i. e. it diverges with m — > , but the velocity structure function has a limit 

lim D ij hc) = Di ((2 + K)5 ij \x\ K - Kx l x j \x\ K - 2 ) (6) 

which is a homogeneous function of x . D\ = 2 2 + K Tr :i ^K(3+^/r((3+K)/2) ^o an d both con- 



stants have dimension 



length 2 
time 



We would like to study the statistical properties of the solutions of Eq. (P in the 
regime of small u, small m (which may be viewed as the inverse of another integral 
scale) and large L. In particular, the universality question for the passive scalar may 
be formulated as inquiring about the existence of the limit of the correlation functions 
< Y\T(t n ,x n ) > in a stationary state of the system when is, m, L^ 1 — > and about the 

n 

independence of such a limit of the shape of the source covariance C . We will show that 
the model possesses an "inertial" range of scales (is/Di) 1 ^ \x\ min(L, m _1 ) where 
these correlators become independent of is, have a limit as v — > and m — > (independent 
of the order), but in general have non-universal (i.e. dependent on the forcing covariance) 
contributions involving positive powers of L. In particular we show that the structure 
functions 

S 2M (x) =<(T(x)-T(0)) 2M >~ 72M(VMr M M (2 ~ K)M (7) 

for \x\ L in the is = limit. The amplitudes are k- and C-dependent and the 
anomalous exponents p<iu depend on k but not on C. We find that p 2 = but 

p 4 = i« + 0(«3) (8) 



for small k. The Holder inequality implies that is a convex function of N . It follows 
that all p 2 M for M = 2,3, ... are strictly positive (p 2 M > (M — l)p4) and that they 
increase with M. Thus structure functions of order four and higher exhibit anomalous 
scaling and have explicit integral scale dependence. While we are able to calculate p 4 
only for small k there is no reason to doubt the generality of the phenomenon for all 
< k < 2. We want to stress that the non-universality is due to the pumping scale, 



not to the dissipation scale, although the derivatives of the field T have short distance 
(non-anomalous) singularities that are easy to analyze ||. 

It is quite well known that in the stationary state the scalar correlations satisfy linear 
PDE's || 0. In the presence of the UV and IR cutoffs v and m, L they have well defined 
representations in terms of the Green functions of the corresponding differential operators 
which we now recall (for more details see || |10|1 ). 

Suppressing the spatial variable, the solution of Eq. ([!]) with the initial condition T 
at t = to takes the form 

TOO = R{t, to) T + /*R(t, s) f{s) ds , (9) 

where R(t,to) is given by the time ordered exponential (t > t ) 

R(t,t ) = rJ> A+v{T) ^ dT . (io) 

Thus, to calculate the correlations of T we need to evaluate expectations of products of 
matrix elements of R(t,to) . We shall use the tensor product notation R(t,to)® N as a 
bookkeeping device for all such products. One then calculates 

< R(t,t ) m > = e -^- to)MN , (11) 

where M. ^ is the differential operator 

N 

Mn = -£(z/A Xn + |r^(0R^) _ ^^( Xn - Xn ,)4,^ . (12) 

n=l n<n' 

The gaussian integral of the time ordered exponentials is calculable due to the inde- 
pendence of v's at different times. The ^T>^(0) term is the contribution of contrac- 
tions within a single R and the last terms come from contractions between different i?'s. 
The former one is an "eddy-diffusivity" contribution: since T> u (0) = ^5^V ll (0) and 
Mi = -{v + £X>"(0)) A , it follows from Eq. ([11]) with N = 1 that, in the absence of the 
sources, the expectation value of the scalar diffuses with the effective diffusion constant 
composed of the molecular diffusivity v and the eddy diffusivity |P"(0) . For small m, 
the eddy diffusivity dominates and the diffusion is driven by the large distance scales 
(recall from Eq. <g) that £>"(0) = O(D m~ K ) ). 

To get hold of the steady state of the scalar, let us first consider the 2-point function. 
From ([5]), we obtain 

< T(tf 2 > = e - {t - to)M2 T(t f 2 + Jds e -(*-^ 2 C L . (13) 

When t - ► — oo, the term with T(i ) disappears due to the positivity of M.2 and we 
obtain for the steady state 

< rp®2 > _ M -1 Cl . (14) 

Due to the translational invariance of C the eddy diffusivity will not contribute to (Q). 
Indeed, note that A4 2 commutes with (three-dimensional) translations and in the action 
on translation-invariant functions of Xx — x 2 = x reduces to 



M 2 = -2uA - V^i^jdidj 



(15) 



Since P y (x) = P IJ (0)— D lJ (x) has an m — * limit given by (|B|), so does the operator A4 2 
in the action on translation-invariant functions and when v — > , it becomes a singular 
elliptic operator Ai s 2 = — D\ ((2 + k) 5 lJ \x\ K — nx l x^ \x\ K ~ 2 ) didj . It is now easy to 
analyze flU]) as the various cutoffs v , m , L are removed using the rotational invariance of 
M. 2 . In the m —>■ and v — > limits (which commute, we could also take m proportional 
to L -1 with no loss), one obtains for the two point function .F2(|x|) =< T(x)T(0) > 

KM = 1,L— - r 2 -(l + OQ) , (16) 

where 72 is a non-universal (i. e. C-dependent) constant and e = |C(0) is the energy 
dissipation rate of the scalar. Note that the non-universal term (a constant) is annihilated 
by .M 2 C . This has to be so if the equation A4 2 ic F 2 = Cl is to be satisfied: the right hand 
side becomes universal in the limit L — > 00 so all non- universal terms in F 2 {r) surviving 
in this limit have to be annihilated by A4 2 c . We shall see this general mechanism limiting 
possible non-universal terms also for the higher point functions. The constant term of 
F 2 drops out from the 2 nd structure function which has a universal L — > 00 limit so that 
the exponent p 2 = 0. The same universal result holds approximately in the whole inertial 
range ?/ < r < min(L, m _1 ) , where the Kolmogorov scale r] = {y jD\^l K . 

Let us now analyze the higher point correlators. Proceeding as with the 2-point 
function, the steady state solution in terms of the operators M.n follows after some 
simple algebra. For the 4-point function one gets 
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< ]J T(x n ) > = F 4 (xi, x 2 , x 3 , x 4 ) + F 4 (xi, x 3 , x 2 , x 4 ) + F 4 (x 1; x 4 , x 2 , x 3 ) (17) 

n=l 

with the single channel function 



F± = Ml 1 ( M 2 X <g> 1 + 1 <g> M 2 Y ) C L ® C L . (18) 

Similarly, the higher stationary state equal-time correlation functions < T® 2M > are 
obtained by symmetrizing the expressions 

M 2 - K }{M 2 - K l 2 ® 1 2 ) • • • (M, 1 ® l 2M - 2 ) C L ® M (19) 

in their arguments (the odd-point functions vanish). Expressions fll8|) and (|l^) are well 
defined for u, m, L~ x nonzero and we need to discuss their limits as these cutoffs are 
removed. 

The main points of this analysis are the following. Acting on translation invariant 
functions as in (|19|), M. N becomes 

N 

M N = A x „ + E ^'( x « - x «') d <% ( 2 °) 

n=l n<n' 



i.e. the eddy diffusivity cancels and (|20[) has an m — > and v — > limit which is a 
singular elliptic operator A4 S § . It can be shown (0) that the Green functions occurring 
in (|I9|) also have limits that are well defined in the UV and render ( |19|) finite for L finite. 
Thus we need to find the leading behavior of dl9|) as L — >• 00 with M. 2 m replaced by 



AA SC 

JV1 2M ■ 



Let us look in more detail at the 4-point function. Recalling that F 2 = AA^Cl , 
it is convenient to view Eq. (|1|) as a differential equation for F 4 that becomes for the 
connected part F 4 = F 4 — F 2 ® F 2 

M4 Fl = £ (F 2 ® F 2 ) , (21) 



where £ is given by the sum in with n = 1,2 and n' = 3,4. By flTBD, the RHS of 
( plf ) has a well defined limit as L — > oo given by 



^^M 2 ^ ( 22 ) 

and is a homogeneous (rotationally-invariant) function of X = xi — X2, Y = X2 — X3 and 
Z = X3 — X4 of degree 2 — k. It is not difficult to write down a solution F£ sc of Eq. (f2~ll) 
for the limiting case with Ai s £ and ( ^2|) on the RHS. One easily checks that 

77CSC e 2 / |v|2(2-k) 1 |r7|2(2~/t)\ e 2 | v 1 2— k I rj I 2-k /oq\ 

^4 = 6(2- K ) 2 (5-k)D 2 '+1^1' ') ~ 9(2 - K )2 D 2 |X| |Z| , (23) 

is such a solution with the use of the decomposition M. s £ = .M| c <8> 1 + 1 <S> .Mf c — £ since 
A4f <S> 1 + 1 <E> -M| c vanishes in the action on F 4 CSC and £ annihilates functions depending 
only on X or only on Z. Thus, we deduce that A4f (F 4 C — F 4 CSC ) — > as L — > oo. By 
scale invariance it is thus reasonable to conclude that the solution for finite but large L 
should differ from the universal scaling form by zero modes of M!f so that 

F c £ L^Y,^ m Fl nm (24) 

0<p4,n<2(2-«) m 

where F£ nm are homogeneous zero modes of Ai s £ of degree 2(2 — k) — p 4i „ and the 
non-universal coefficients 7„ m depend on the source covariance C . 

In fact using spectral analysis of M.± 0, (|24| ) can be made rigorous (possibly with 
logarithmic corrections in L for special values of k) . Similar analysis can be repeated for 
iV-point correlators: non-universal L-dependent terms proportional to homogeneous zero 
modes of M. s § can be present in the large L asymptotics. We thus face the problem of 
finding such zero modes, of determining whether they are present in the iV-point function 
of T and finally of finding whether they contribute to the structure function Sjy While 
such an analysis still eludes us for general values of k, we will now show that at least for 
small k the zero modes are present and dominate the structure functions. 



The stationary state correlation functions of the scalar, given in terms of (|19f) , are 
in general non-gaussian, but they become gaussian as k — > 0. To see this, note that in 
this limit D %: > = 2 Di 5 iJ (having finite Di requires the vanishing of D Q as K — > in 
order to renormalize the ultraviolet divergence in @; Dq will never show up below.) We 
immediately obtain for the k = operators : 

Mf tQ = 2AV X1 • V X2 = -2D, A x , (25) 

l<n<n'<4 

= -2D 1 (A x + A Y + A z - V x • V Y - V Y ■ V z ) (26) 

in the difference variables X, Y, Z and using the subscript to refer to k — 0. Some 
straightforward algebra shows that the expression (|19D reduces (when symmetrized) 
to the standard gaussian expression of sums of product of 2-point functions. 

Our strategy is now the following. We shall find the homogeneous zero modes of the 
operator Aif in perturbation expansion in powers of k. Eq. (^) implies that (for m — 0) 

5^(x) = 2Di(<J« + k i? J '(x)) + 0(k 2 ) with 

iTJ'(x) = 5 ij (^ + In |x| ) - ixV" |xp 2 (27) 



Hence, to the first order in k, M.*f = A4^ + 2kD\ V4 , with 

V 4 = -ffi(X)d X id X i - ir{Y)i)y.<) y . - R ij (Z)d Z id Z i 

- (R ij {X + Y) — R^iX) - R ij {Y)) d xi d Y o 

- (R ij {Y + Z) - R ij (Y) - R ij (Z)) dyi d zi 

- (R lj (X + Y + Z) - R ij (X + Y) - R lj (Y + Z) + R lj (Y)) d x > d zi (28) 

Since A4| c commutes with 3-dimensional translations and rotations and with the 
permutations of four points we shall search for its zero modes respecting these symmetries. 
Consider first the unperturbed zero modes. The symmetric zero modes of the lowest 
homogeneity of A^| c occur in degree zero (constants) and in degree 4. The latter, when 
expressed in terms of point-differences, have the form 

a^(x n -x n /) 4 + o^(x n -x m ) 2 (x„-x m /) 2 + cj^ (x„-x n /) 2 (x m -x m /) 2 

{n,n'} {{n,m},{n,m'}} {{n,n'},{m,m'}} 

= aF 1 + bF 2 + cF 3 , (29) 

where the pairs {n,n'} and {m,m'} are assumed different, as well as the pairs {n,m} 
and {n, m'} and where 10a + 146 + 3c = . 

The constant survives as the eigenvalue of Ai s £ for Thus we need to calculate in 

degenerate perturbation theory how the fourth degree zero modes change with k. For this 
we write Mf >0 = -2£>i ( + + A^), where X = X, Y = ^(Y + ^X+^Z) , Z = 
Z . Denoting R = (X 2 + Y 2 + Z 2 ) 1/2 = (J £ (x„ - x n ,) 2 ) 1/2 , we obtain 

{n,n'} 

Ml, = - 2 -§-d R R 8 d R - ^$ (30) 

where $ is the Laplacian on the sphere S 8 in the space of (X, Y, Z) . Now pick two lin- 
early independent zero modes R /» , i — 1, 2 , of the form ( ]2"9"| ) and look for a homogeneous 
zero mode 

R A+K \a l f l + a 2 f 2 + nh) (31) 

with a homogeneous degree zero function f 3 orthogonal to /i j2 in L 2 (S 8 ) . We obtain 
in the linear order in k 

MT fi (\RUnR(a 1 f 1 + a 2 f 2 ) + R A / 3 ) + 2D, V 4 R 4 (aj, + a 2 f 2 ) = (32) 
or, using the form ( p0|) of M s £ , 

- 15 X(a 1 f 1 + a 2 f 2) -44/3 - $/ 3 + ^^^(oiA + aa/a) = . (33) 
Upon taking the L 2 (S 8 ) scalar products with /i )2 , ^3 drops out resulting in the relation 

-15A E(/*.^^/i)ai = ° ( 34 ) 

3=1,2 ./ 1.2 

Hence A has to solve the equation 

det ((f it ^V.R'fj) - A(/<, £)) = . (35) 



For the explicit calculation we took R A f\ = 3Fi — 10F 3 and R 4 f 2 = —7F\ + 5F 2 . The 
integrals over the 8-dimensional sphere are most conveniently done by using homogeneity 
to transform them to gaussian integrals over R 9 . The integration is straightforward with 

MapleQ and the matrix in fl3~5f) becomes proportional to ( ig^ — 20~— ^20A ) 

the eigenvalues Ai = —14/5 and A2 = —1. The corresponding eigenf unctions fl3T|) are 
given by F 4 C 1 = R-^+^) K +°^ 2 )(F l - 2F 2 + 6F 3 + 0(k 2 )) and F 4 C 2 = R~ K+0 ^(-7F 1 + 
5F 2 + O(k)) . For large L the connected 4-point function takes the form 

< f[ AO > C = 74,o L*~ 2k + 74,1 (L/R)^ 5+ °^ R- 2K (Ft - 2F 2 + 6F 3 + O(k)) 

+ 74,2 {L/R)~ K+0 ^ R~ 2K (-7F l + 5F 2 + 0{k)) 
4 

+ < J] T(x n ) > csc + 0((L/R)- 2+0 ^) (36) 



71=1 



n=l 



uniformly in small k. Since the connected correlation vanishes and < n T(x n ) > 

n=l 

2 2 

reduces to 36 q D h (3-Fi — IOF3) for k = 0, we infer that 7^0 = , 74,1 = 2 i6D' 2 + 
and 7 4j 2 = 5i Q D 'i + C(^) • The result (0) for = 2 follows with 74 = 3^ + O(k) and p 4 

given by (^|) since only the F/j^-term gives non-zero contribution to the structure function. 
The zero mode, with A = — 1, is actually obtained from a zero mode of .M3 by extending 
it to a function of four Xj's by symmetrizing. This is a general feature: zero modes of 
A4n give rise to zero modes of A4 2 m for 2M > N. These however do not contribute to 
the structure functions S 2 m- The only zero mode of A4 2 m that contributes, is the unique 
one not coming from the lower dimensional .Mat's, namely the one that at k = is gotten 
from the monomial nf=i( a; 2j-i — x 2 i) 2 by symmetrizing and subtracting partial traces. It 
gives rise to the dominant contribution to S 2 m which has to be present by the Holder 
inequality. 

The asymptotic behavior of the scalar correlation functions encodes a subtle informa- 
tion about the behavior of the Green functions of the singular multibody operators M. n 
with continuous spectrum. The reduction of its study to that of discrete spectrum oper- 
ators given by .Mat's acting on homogeneous functions should be thought of as realizing 
a renormalization group type approach to the model, with the homogeneous zero modes 
of .Mjv playing the role of relevant interactions. This may be the most important hint 
from the above exact solution for the anomalous scaling of the passive scalar. 
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